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Abstract. Strange modes have been found in the radial 
spectrum of many luminous stars, such as PAGB stars. 
The strange modes are characterized by having small am- 
plitudes in the interior of the envelope, and egregious pe- 
riods and growth-rates. It has been common belief that 
the strange modes are a result of strong nonadiabaticity. 
Recently, and perhaps surprisingly, such modes have also 
been found in classical Cepheid models, even though these 
are weakly nonadiabatic stars. Here we show that in fact 
there is nothing strange about these modes and that they 
must exist even in the adiabatic limit. They are essentially 
acoustic surface modes. 

By means of a simple change of variables and without 
approximation, the adiabatic linear pulsation equation for 
the radial displacement is reduced to a Schrodinger like 
equation in which the radial coordinate is the local sound 
traversal time. In this formulation, the narrow hydrogen 
partial ionization region is seen to act as a potential bar- 
rier, separating the star into two regions. Modes can be 
trapped either in the inner or in the surface region. Cou- 
pling through the barrier gives rise to resonances between 
the inner and surface regions. The strange modes are those 
in which the ratio of inner to surface amplitude is at a 
minimum. The potential problem formulation shows that 
strange modes exist in the adiabatic limit. As a function of 
the stellar parameters the relative location of the barrier 
changes, and this gives rise to the phenomenon of avoided 
level crossings along a sequence of models. 

The appearance of strange modes and the associated 
level crossings can be exhibited with an analytically solv- 
able toy model when the potential barrier is approximated 
by a delta function. In the full nonadiabatic models the 
same trapping mechanism remains responsible for the ap- 
pearance of strange modes. The unusual growth-rates are 
seen to also be a consequence of the relative minimum of 
the inner amplitude for these modes. Again the behav- 
ior of the nonadiabatic modes can be well mimicked by a 
simple analytical toy model. 
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The strange modes can be linearly unstable to the left 
of the fundamental and first overtone blue edges. Hydro- 
dynamical calculations show that the strange limit cycle 
pulsations (a) are extremely superficial as the linear eigen- 
vectors already indicate, in fact they have negligible am- 
plitudes interior to the partial hydrogen ionization front, 
and (b) the pulsations have surface radial velocities in the 
0.1 - 1.0 km/s range, but extremely small photosphcric ve- 
locities, and luminosity variations in the milli-magnitude 
range. These modes are therefore expected to be difficult 
to observe. 

Key words: stars: variables: Cepheids - stars: oscillations 



1. Introduction 

This work sets out to demystify the so-called strange vi- 
brational modes that are regularly encountered in stars 
with high luminosity to mass ratio, such as PAGB stars 
(Cox et al. 1980, Saio et al. 1984, Aikawa 1985, Gautschy 
& Glatzel 1990, Zalewski 1992, 1994, Cox & Guzik 1996, 
Papaloizou et al. 1997), but more recently also in Cepheids 
(Buchler 1997) where they were not expected because of 
the weakly nonadiabatic nature of the pulsations in these 
stars. We begin in §2 by rephrasing the description of adia- 
batic pulsations in the general form of a quantum mechan- 
ical potential well problem. The natural follow-up is to 
focus on the resulting potential. This we do for a Cepheid 
and find that the sharpness of the hydrogen ionization 
front produces an immensely high barrier in the potential, 
effectively separating the star into two distinct regions. 
The barrier is ultimately caused by the rapid variation 
of the local sound speed, most by the adiabatic index. In 
§3, for didactic purposes, we present a toy problem that 
approximates the Cepheid, yet contains all the essential 
features of strange modes, and that can be solved ana- 
lytically. We thereby verify that the first appearance of a 
strange mode occurs in concert with a node crossing the 
potential barrier into the narrow outer envelope of the star 
between the surface and the hydrogen partial ionization 
region (HPIR). 
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In §4 we look in detail at a specific Cepheid model, us- 
ing the transformations inspired by the potential well for- 
mulation set out in §2. An entire sequence of Cepheid mod- 
els is then examined in §5, and the signature of strange 
modes is revealed unequivocally. It will become apparent 
that the gradual change in the model parameter, which 
affects the extent of the region between the HPIR and the 
stellar surface, modifies the resonance conditions between 
the inner and outer regions and is responsible for the ob- 
served systematic changes along the sequence. In §6 we 
extend the adiabatic toy model of §3 to the nonadiabatic 
case by allowing a complex potential barrier. Although we 
show this model to be highly approximate, it does repro- 
(lucx; the behavior of the periods and of the growth-rates 
in Cepheid sequences quite accurately. 

The pulsating stellar envelope, both physically and 
mathematically, bears a strong resemblance to a (quar- 
ter wave length) flared wind instrument such as a horn 
or trumpet. In §7 we explore this analogy and design an 
equivalent musical wind instrument with the characteris- 
tics of the Cepheid spectrum. 

The strange modes can be unstable outside the usual 
Cepheid instability strip and can thus give rise to nonlin- 
ear pulsations of their own. In order to shed some light on 
the appearance of such pulsations and answer the question 
of their observability we therefore compute the nonlinear 
hydro dynamical behavior of several such "Cepheid" mod- 
els in §8. 

2. The Cepheid Potential 

In the standard equilibrium radiative diffusion approxi- 
mation adiabatic radial pulsations arc described by the 
familiar equation for the radial displacement eigenfunc- 
tions Sr: 
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where the symbols have their usual meanings (e.g. Cox 
1980). 

We can recast the description of the adiabatic pulsa- 
tion problem in terms of a new function $ by invoking the 

transformation 
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Fig. 1. Effective potential for the adiabatic pulsations; top: 
full potential (solid line) and gravitational contribution (dot- 
ted line); the horizontal lines denote the lowest eigenvalues w^; 
middle: scale transformation r{z); bottom: blow-up of the re- 
gion near the hydrogen partial ionization region; Stellar surface 
is on the right at « = 1. 



$ = (riP)2 rSr = (p)2 CsrSr , 



(2) 



where the second equality simply involves the adiabatic 
sound speed: = TiP/p. Thus 



(3) 



where the assumed exp{—iujt) time dependence has al- 
lowed us to replace the temporal derivative by the eigen- 
value w^, and the function W{r) is given by 



W{r) = 



(4) 



The first term of W{r) embodies the effects of the spa- 
tially varying inertia (density) and sound speed on and 
the form it takes is a result of applying Eq. (2) to Eq. (1). 
Anywhere the equilibrium varies rapidly. W{r) is domi- 
nated by this first term. The second term of W{r) results 
from the variation of gravitational forces within the hy- 
drostatic equilibrium and will only be important at small 
radius and for modes with small wave number (these have 
local wavelengths comparable to the scale for which the 
gravitation varies). Both terms play a role in determining 
the global frequency spectrum of radial modes. 
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We note in passing that when a plane-wave function 
$(r) = ex.p{ikr) is inserted, Eq. (3) then represents the 
local dispersion relation and can be used to construct a 
propagation diagram (Unno et al. 1989). 

We can now go one step further and convert Eq. (3) 
into a Schrodinger like equation by eliminating the sound 
speed that modifies the spatial derivative in Eq. (3); we do 
this by transforming from r to a new coordinate z, chosen 
to increase inward: 

dz = —dr/cg. (5) 

The quantity z then represents the acoustic depth (or 
sound traversal time) and is conveniently measured in 
sound days here {r{z) appears in Fig. 1, center panel). 
The newly phrased problem involves a concomitant trans- 
formation of the eigenfunctions according to: 

^ = {pCs)^rSr = Cs (6) 
which converts Eq. (1) into 

-^^ + V{z)^ = ujH, (7) 
where the potential is given by 

Note that the potential V in Eq. 7 differs from W. 

Eq. (7) is clearly the simplest and most compact form 
that the problem of adiabatic stellar pulsations can be ex- 
pressed in. The spatial gradients of the equilibrium quan- 
tities p and Cg, as well as sphericity effects have been in- 
corporated into the modified eigenfunctions \1/ and into 
the transformed spatial coordinate z. The effective sound 
speed is therefore unity and wc can set lu = k, the wave- 
vector. The boundary conditions associated with the orig- 
inal problem Eq. (1) are readily transcribed to the new 
formulation of Eq. (7). 

Eq. (7) has been derived for the radial displacement 
eigenvector, or equivalently the velocity. In the litera- 
ture one also finds equations similar to Eq. (3) or (7), 
but for the pressure eigenfunction 6p (e.g. Deubner & 
Gough 1984). We stress that the derivation for such a 
Schrodinger equation for Sp is only possible if approxima- 
tions are made. The exact equation for dp would involve 
third order spatial derivatives, unless the variations of the 
gravitational potential are neglected in the linearization 
(Cowling approximation). 

We now examine the potential of a specific Cepheid 
model, choosing the parameters L = SOOOLq, M = SM©, 
Te//=5520K, Z = 0.02, X = 0.70. We use the latest 
OPAL opacities (Iglesias & Rogers 1996) with Alexander- 
Ferguson (1994) molecular opacities. Convection is ig- 
nored. 



The top panel of Fig. 1 displays V{z) (heavy solid 
line) and its gravitational component (dotted line) for this 
Cepheid model. The gravitation contributes a small, but 
significant depth in the interior. In this plot the coordi- 
nate z ranges over the entire pulsating envelope. The scale 
of z is of the same order as the period of the gravest pul- 
sational modes (cf. virial theorem). In order to reveal the 
overall structure of the potential well the vertical scale is 
highly magnified. The horizontal lines indicate the relative 
placement of the adiabatic pulsation energies. The barrier 
at 2 = 0"? 75 is so high as to be off-scale and so narrow as 
to be featureless with respect to the scales that are appro- 
priate for the pulsations. Therefore in the bottom panel 
of Fig. 1 we exhibit the shape of the potential V{z), re- 
stricting our perspective here to a very narrow region in 
the vicinity of the HPIR. 
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Fig. 2a. Periods of the barrier problem, with 14=10; Zo is the 
distance of the barrier from the inner (nodal) boundary; 1-Zo 
is thus the width of the outer region. 
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Fig. 2b. Period ratios of the barrier problem, with 14=10; 1-Zo 
is the width of the outer region. 
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Fig. 3. Lowest four eigenfunctions as a function of barrier lo- 
cation Zo, with 14=10; 1-Zo is the width of the outer region. 



The middle panel of Fig. 1 plots the relation between 
the radial coordinate r and the acoustic depth z. In the 
new z coordinate the inner regions with their high sound 
speed now appear highly compressed. We have adopted 
zero acoustic depth {z = 0) at the stellar surface. 



Of the three quantities that appear in the second 
derivative of Eq. 8, it is the sound speed whose variations 
are most significant. Recall that oc and F expe- 
riences rapid dips (e.g. Cox & Giuli 1968) in ionization 
zones that then produce the large potential fluctuations. 
The centers of the three ionization stages, H-H+, He-He"*", 
He"'"-He"'""'" are indicated by arrows in Fig. 1. Of these it 
is the hydrogen that is by far the narrowest, leading to a 
potential barrier. 

3. An Idealized Potential Barrier Toy Problem 

For a potential like that shown in Fig. 1 the eigensolu- 
tions have many properties that can be demonstrated in 
a more transparent way by looking at a toy problem. This 
toy problem can be solved analytically, yet it displays the 
essential features of the modal spectrum of the star. We 
replace the actual potential with an infinitesimally narrow 
barrier (a delta function of strength Vo) which is located at 
Zo (0 < < 1) and the potential simply zero elsewhere. 
Thus the Schrodinger equation reduces to 



-y" + Vo k6{z - Zo}y = iJ^y. 



(9) 



where k5{z — Zo) models the potential shape and Vo is 
a dimcnsionlcss scale factor. Wc choose the boundary 
conditions appropriate for the adiabatic pulsations of a 
Cepheid, those for a quarter wave-length tube that is 
closed at the inner boundary and open at the outer (pure 
reflection), viz. : 



y(0) = 2/'(l) = . 

The eigenvectors are given by 

y{z) = a sin kz 

= cos fc(l — z) 



Z < Zo 
Zo< Z 



(10) 



(11) 



with normalization y(l) = 1 (more generally a{k) = 
y'{0)/ky{l)). With the appropriate matching conditions 



y'{z+) = y'{z-) + kVoy{zo) 



15/16 one obtains the characteristic equation for k 
cos A; -h 14 sin kzo cos k{l — Zo) = . 



The inner amplitude a is given by 
a{k) = cosA;(l — 2:o)/sinfc^;o- 



(12a) 
(126) 



(13) 



(14) 



We note several interesting features. First, Eq. 13 
shows that for large Vo (a very high barrier) one obtains 
cos A;(l — Zo) sinfczo « 0. Thus for Zo ^ (vanishing inner 
region) the spectrum is given by cosfc = 0, with periods 
Pn = 4/(2n + 1), i.e. the one of an open tube, whereas 
for Zo ^ 1 (vanishing outer region) we have essentially a 
closed tube and its spectrum sin fc = with P„ rs 2/ (n+1). 
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4. Eigenvectors j/n of the potential barrier toy problem 
Zo = 0.96 and with Vo=W; the strange mode is labelled 



This asymptotic behavior in Zo is nicely seen in Fig. 2a 
that shows the period spectrum as a function of 1 Zq for a 
value 14=10. (The reason for plotting versus I-Zq instead 
of Zo is that we eventually compare to fig. 10 in which 
the outer acoustic distance increases with Pq.) Following 
common practice we assume that the modes are ordered 
by decreasing periods P„. 

A null inner amplitude can occur for specific combina- 
tions of Vo and Zo, i.e. nongenerically, as Eq. (14) shows. 
We call these modes perfect surface modes. It will be seen 
later in the application to stars that the surface modes are 
in fact the strange modes. For arbitrary values of Vo and 
Zo we refer more loosely to the mode whose amplitude is 
at a minimum as the strange mode. 

The outer boundary condition precludes a zero outer 
amplitude, i.e. infinite values of a, but a can get very large. 
We refer to the corresponding modes as inner modes; later 



we shall see that in stars the inner modes are the regular 
modes. 

Finally, in between the inner and the outer modes lie 
modes that have a uniform amplitude throughout (a=l). 
Such perfectly matched modes can occur accidentally 
when a node happens exactly at Zq (cf. Eq. 12b). For most 
values of Vo and Zo one can generally identify modes that 
are close to being matched. 

It is amusing to notice in Figs. 2aand 2b how, through 
the intermediary of trapped modes, the spectrum gradu- 
ally changes from that of a closed to a open tube as I-Zq 
varies between and 1. The reason for this, as we have 
seen, is that a very strong barrier (large Vq) acts almost 
like an additional interior boundary condition. Below the 
barrier, the modes resemble those of a closed pipe of length 
Zo, while above the barrier the modes resemble those of 
a open pipe of length l-Zo- Except for the rare situation 
of matched modes, either the inner or the outer ampli- 
tude dominates, i.e. the modes are either inner (regular) 
or outer (strange) modes. 

The appearance and behavior of the modes is further, 
and more globally illustrated in Fig. 3 which shows the 
lowest three eigenvectors as a function of z. We recall that 
the closed end is on the left {z = 0) and the open end 
on the right (z = 1). In the right margin we show the 
location of the barrier Zq. This figure should be looked at 
in conjunction with Fig. 2a which shows the periods P„ 
as a function of Zq- 

Consider first the fundamental mode displayed in the 
leftmost column of Fig. 3. As expected, for small Zq the 
outer well dominates and the eigenfunction is essentially 
that of an open, quarter-wavelength tube. As the barrier 
moves to the right, the effective length of the tube shrinks, 
causing a decrease in the period as seen in Fig. 2a. On the 
other hand, starting from the bottom, we see that the 
eigenfunction is now essentially that of a closed half wave- 
length tube with half the period of the Zq = case;. For 
Zo=2/3 one can fit exactly 3/4 of a wavelength into the 
well, which corresponds to a perfect match condition. At 
the shown values of Zo=10/16 and 11/16 the yo eigen- 
vectors axe indeed very close to those of the first over- 
tone of an open tube, but being the fundamental mode 
they avoid acquiring a node. This is also the place where 
a near degeneracy occurs as Fig. 2a indicates. However, 
because the Sturm-Liouville nature of the problem pre- 
cludes degeneracy, there must be an avoided level cross- 
ing. (Strictly speaking, the problem is not Sturm-Liouville 
because of the 5 function potential, but seems to preserve 
the features thereof in this singular limit). If we look at the 
second column which represents the first overtone in the 
well we notice that the system avoids degeneracy through 
a switchover between the fundamental and the first over- 
tone, indicated in Fig. 3 by the X linking columns 1 and 2. 
There is of course nothing special about the fundamental 
and first overtone, and we see similar avoided crossings in 
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Fig. 4 and corresponding switchovers in the eigenvectors 
in Fig. 3 that are also indicated by crosses. 

The lines with a positive slope in Fig. 2a correspond 
to modes with predominant amplitudes in the outer well, 
i.e. with open tube characteristics, whereas the those with 
negative slope have predominant amplitudes in the inner 
well, i.e. are essentially closed tube modes. It is thus the 
overlap of the open and closed tube spectra that makes 
up the pattern of Fig. 2a. 

The shapes of the fundamental eigenvectors in Fig. 3 
show that changing Zo is essentially equivalent to chang- 
ing the length of the tube, and explain why the periods 
increase/decrease linearly with Zo- 

The periods P„ change with Zo^ but the ratios P„o arc 
almost everywhere independent of Zo and thus straight 
lines. The reason for this near constancy is again that 
changing Zo changes the equivalent length of the tube. We 
show the period ratios in Fig. 2b so that later they can be 
compared to those of the Cepheid model sequence. 

Why do the modes alternate between inner and surface 
modes? Consider again Fig. 3, column 1 starting at the 
bottom and going up. Near the open end the eigenvector 
is forced to bend down, at a constant rate because the 
curvature of y is given by —k'^y. Therefore the slope of 
the eigenvector in the outer region gets increasingly steep 
as the barrier location moves inward. Since there must be a 
finite, fixed change in slope at Zq (Eq. 12b) the slope of the 
eigenvector of the inner region is increasingly flattened, 
thus suppressing the inner amplitude (which has to go 
to zero at Zo=0)- The overtone modes follow the same 
scenario. 

We now examine more closely the case with parameter 
values Zo = 0.96 and Vq = 10 that approximately mimics 
the Cepheid situation. Eq. (13) is easily solved and the the 
lowest 14 modes eigenvectors yn are displayed in Fig. 4 vs. 
acoustic depth z. Referring back we see that this puts us 
on the extreme left of Fig. 2. We note in particular that 
the strange (surface) mode is expected to appear only for 
large mode number. Indeed, Fig. 4 shows that here it is 
mode Oil that is strange (and mode 012 that is nearly 
matched). Note that for mode 012 a node occurs very 
close to Zo which then causes no jump in the slope (Eq. 13). 
Modes higher than 012 therefore have at least one node 
in Zo< z< 1. 

One can again look at the behavior of the modes in 
Fig. 4 in terms of matching of the eigenvectors through 
the potential barrier. Here now the curvature increases 
with modal index (wave number k). For the lowest modes 
the curvatures and bending down of the eigenvectors are 
small and the modes arc inner (regular) modes. As mode 
09 is approached the matching has already started to se- 
riously suppress the inner amplitude which reaches a min- 
imum for Oil which is the surface (strange) mode. Note 
that both OlO and 012 have approximately the same 
inner and outer amplitudes and are hybrid modes. For 
mode 012 furthermore the node occurs very close to the 



barrier location Zq and 012 is thus an almost perfectly 
matched mode. The subsequent modes become again reg- 
ular modes, and through the same mechanism a new sur- 
face mode appears. This whole scenario repeats itself with 
a relative recurrence period of w 1/(1 — Zq) in the modal 
index (cf. Eq. 14). 

As a byproduct of the previous paragraph we note that 
the minimal amplitude and the near match occur very 
close to each other, and that both phenomena are the 
consequence of a near resonance between the inner and 
outer potential wells. 

The connection with the usual stellar modes will be- 
come more apparent when we consider the behavior of the 
Svn rather than the 'i'n in the next section. 

Thus in the Cepheid case the radial vibrational modes 
that are confined to the surface regions have strange prop- 
erties. Under cross-over conditions both modes are an ad- 
mixture of normal and strange modes. This is of course not 
all that different from what happens for nonradial Main 
Sequence stars where different mode types, g and p modes, 
are involved in the level crossings (e.g. Unno et al. 1989). 

4. Details of a Cepheid Model 

Before discussing the strange modes in sequences of 
nonadiabatic Cepheid models where they were first no- 
ticed (Buchler 1997) we first return to the Cepheid model 
of §2 in order to show how the toy problem of §3 serves as 
a perfect paradigm for explaining the nature of the strange 
mode. 

We have to be a little bit careful with the outer bound- 
ary condition because of the almost isothermal nature of 
the outer region. In order to avoid problems we define 
the surface as the point where the density is small, but 
nonzero, or more conveniently where the pressure is equal 
to some small value , say two times the radiation pres- 
sure (aT^/3). We impose this as a constant outer pres- 
sure in the subsequent dynamics. (The exact value of this 
pressure is not critical). 

4-1- The Adiabatic Cepheid Model 

A great deal can be learned by examining the linear pul- 
sations in the adiabatic limit, i.e. we ignore the entropy 
variations during the pulsation. 

The adiabatic displacement eigenfunctions for the 
first fourteen pulsational modes are shown in Fig. 5, start- 
ing with the fundamental mode (F). With increasing k a 
node moves outward toward Zo and, as overtone mode Oil 
is approached, the node is about to reach Zo- The similar- 
ity with the analytical toy-model of the previous section 
is striking. As in the toy problem, the inner amplitude de- 
creases to a minimum, - the property that has been used 
to identify a strange mode. There is no doubt that in an 
actual Cepheid model the occurrence of the strange mode 
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Fig. 5. Adiabatic eigenvectors of the scaled displacement Fig. 6. Adiabatic eigenvectors 5rn/r vs. acoustic depth for the 

vs. acoustic depth for the Cepheid Model; the strange mode is Cepheid Model; note the different scales on the left and on the 



labelled by <S. 



right column. The strange mode is labelled by S. 



is also the result of the confinement of the mode to the 
outer cavity of the potential well. 

There are some small differences with the toy prob- 
lem, for example in the detailed behavior both near z =1 
and through the potential barrier. This is not astonishing 
considering that in the toy model we ignored the the po- 
tential dip near z = Zo s& well as the broad, but relatively 
low potential barrier near z = 1 (The effect of the latter 
is seen in the exponential behavior of the low modes near 
the surface. 

It is interesting of course to compare the scaled eigen- 
functions to the radial eigenf unctions i5r„ which are 
shown in Fig. 6. The scale factor {r^pcs)~^ in Eq. 6 
greatly magnifies the surface region and the strange mode 
now appears even more as a surface trapped mode. 

While the radial displacement eigenfunction (5r„ has 
perhaps a more physical meaning, the representation 



provides an optimal scaling of the problem, and is perhaps 
more basic from a mathematical point of view. It shows 
in particular that even for a strange mode the inner re- 
gions play an important role even though in 5rn all modes 
appear to be surface modes. 

We conclude that a strange adiabatic mode has its ori- 
gin in the potential barrier that is created by the ionization 
front of hydrogen in the outer envelope. It first appears 
when, with increasing mode number, an actual node is 
forced to occur near the surface (outside the HPIR). This 
then pushes down the inner amplitude of the mode to 
a very low value. The new, interesting point is that the 
strange mode already appears in the adiabatic limit. 
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mode labelled by S; note the different scales on the right and 
left. 



4-2. The Nonadiabatic Model 

We now turn to the correct, nonadiabatic linearization 
of the Cepheid model. In Fig. 7 we display the moduli of 
the nonadiabatic displacement eigenvectors related to 
the radial eigenvectors through the scaling of Eq. (6), and 
normalized to 1 at the surface. Juxtaposed we show in 
Fig. 8 the corresponding radial displacement eigenvectors 
(5r„/r. The nonadiabaticity is weak for Cepheids (k;„ -C w„ 
for the low modes) and therefore the properties of the 
nonadiabatic eigenvectors to a large extent mimic those 
of the adiabatic problem. (We recall that in our notation 
the eigenvalues for an assumed exp(tT„f) dependence are 
denoted by (t„ = iiVn + 



4-3. Transition From Adiabatic to Nonadiabatic Model. 

The nonadiabatic spectrum is noticeably different from 
the adiabatic one as Table 1 shows, in particular it is 
mode 08 which is now the strange mode (denoted by the 
symbol S). It is naturally of interest to explore how the 
transition from adiabatic to nonadiabatic occurs. To that 
effect we show in Fig. 9 the period ratios as a function of 
the logarithm of the nonadiabaticity parameter e (We use 
the same definition as Zalcwski 1992; viz. in the dynam- 
ical equation e multiplies the entropy contribution to the 
pressure variation in the mechanical equation, e =0 thus 
being the adiabatic case and e=l being the full nonadia- 
batic one). 

Note that it is indeed for the strange modes that the 
period ratio changes most rapidly in Fig. 9. The rea- 
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Table 1. Comparison of exact and adiabatic periods 
Prom left to right: mode, period, growth-rate, moment of 
inertia, adiabatic period; the strange mode is indicated 
with S. 
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son for this extreme sensitivity is also clear. The reso- 
nance condition is very sensitive to the outer region. But 
there the Cephcid envelope is highly nonadiabatic (e.g. 
Cox 1980) even though, overall, the Cepheid pulsation is 
weakly nonadiabatic, as already pointed out. It is there- 
fore expected that even a small nonadiabaticity can have 
the dramatic effect of gradually shifting the strange mode 
from Oil down to 08 as Fig. 9 indicates. 

5. Cepheid Sequence 

We now turn to a Cepheid model sequence with M=^Mq, 
i=3OOOi0, X=0.7, Z=0.02, in which the variable param- 
eter is the fundamental period Pq , or equivalently the Teff 
ranging here from 8600 K to 4600 K. We have chosen 
this somewhat overluminous sequence (on the Chiosi M- 
L relation) simply because the strange mode appear as a 
vibrational overtone with lower modal number. 

5.1. Modal spectrum 

We start again by examining the adiabatic models, and 
show in Fig. 10 the behavior of the adiabatic period ra- 
tios P„o = Pu/Pq for the model sequence. The adiabatic 
modes are seen to display avoided level crossings. Here we 
can prove that the levels cannot cross. Indeed, because 
of the Sturm-Liouville nature of the adiabatic eigenvalue 
problem (e.g. Ledoux & Walraven 1958) the modes have 
distinct frequencies - no degeneracy and no crossings. 

Guided by the toy problem (Figure 3) we understand 
that the level crossings are a consequence of the relative 
shift in the location of the potential barrier, and thus 
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Fig. 10. Cepheid model sequence; Period Ratios P„o of the 

adiabatic spectrum based on the first fifteen overtone modes 
of the spectrum (k=l, 15); only ratios with fc > 3 are shown. 

in the resonance condition. Indeed with increasing period 
(or decreasing Te//) the HPIR moves inward, causing the 
strange mode to move to lower modal number. As we go 
from high to low T^ff (low to high Po) the HPIR moves in- 
ward, and concomitantly the outer potential well increases 
in width relative to the inner one. In the toy problem this 
corresponds to decreasing z^. We therefore anticipate a 
behavior comparable to that of the toy problem displayed 
in Fig. 3. 

Turning now to the full noimdiahatic models we plot in 
Fig. 11, on top, the period ratios which also show avoided 
crossings However, they are shifted with respect to the 
nonadiabatic picture and are more dramatic. The mode la- 
belled 15 on the left successively trades places with modes 
14, 13, 12, etc and ends up as mode 7 on the right side of 
the figure. 



10 J. R. Buchler@l, P. A. Yecko@l, Z. Kollath@l @2: The Nature of Strange Modes in Classical Variable Stars 




0.4 
0.2 


-0.2 
-0.4 
-0.6 
0.4 
0.2 


-0.2 
-0.4 



-0.6 - 




-1 1 1 r- 



_i I I L_ 



-1 1 1 r- 




10 

Po (days) 

Fig. 11. Nonadiabatic Cepheid model sequence, top: Period 
ratios, P„o,k=l,14; middle: relative growth-rates, fe=0,6; bot- 
tom: relative growth-rates, k =7, 12; individual models are 
indicated by points in the bottom graph. 

The middle and bottom graphs in Fig. 11 display 
the relative growth-rates rjn = 2KnPn, for modes k=0 
through 6, and k=7 through 14, respectively. Thus the 
first overtone in this sequence is unstable for T^ff < 6000 K 
(Pi =4^7) and the fundamental for T^ffK 5800 K (Pq = 
7^7). 

In the bottom graph that shows the higher overtones 
we witness some interesting behavior. The lines connect 
the points ordered by decreasing P„o which is different 
from the way the eye connects them. There is an unstable 
modal feature throughout the whole period range, but it is 
not always the same mode. Furthermore whenever a level 
crossing occurs in the top graph, there is a corresponding 
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Fig. 12. Work integrands vs. acoustic depth for Cepheid 
Model; the strange mode is labelled by S. 



switch-over in the growth-rates of the two corresponding 
modes. The switchover is relatively slow for modes 7 to 
8 and 8 to 9, more rapid for modes 9 to 10 and 10 to 11 
and very rapid for the higher modes. (The reason is simply 
that the position of the HPIR moves more rapidly with Pg 
for low Pq.) In both the period ratios and in the growth- 
rates it is as if some egregious modal property hops from 
mode 14 to mode 7. 

A more careful inspection of Fig. 11 uncovers among 
the higher modes the appearance of a second strange mode 
in the period ratios. Some secondary structure is even 
more clearly visible in the growth-rates. There are local 
maxima for modes 13 and 14 at 10^5 and 8'}5 days, re- 
spectively, which are due the appearance of higher order 
strange modes - by which it is meant that a second node 
crosses the HPIR. 
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In a recent paper Glasner & Buchler (1993) already 
noted that for RR Lyrae and Cepheid models the stability 
of the vibrational modes does not increase monotonically 
with mode number. Instead, there is a strong excursion 
back toward instability around the eighth or ninth over- 
tone. At the time they had no physical explanation for the 
nonmonotone behavior, but the present results show that 
this is related to the recurrent appearance with increasing 
modal number of strange modes. 

5.2. Work-integrands 

Fig. 11 shows that the strange mode has a very differ- 
ent growth-rate from its neighbors and it is interesting 
to inquire why. In Fig. 12 we therefore examine the work 
integrands, i.e. the work (pdv) done per cycle, or the in- 
tegrand of the expression for the relative growth-rate 



-. 2— = '^^ Im [ SpnSv* 4:TTr'^ pCs dz 

i^n i^iln J 

= -^y- / ApAvsm{(l>p - (j)y) Anr"^ pCs 



dz 



(15) 



where the 4>p is the phase of 6pn and Ap„ its modulus, and 
similarly for the specific volume eigenvector 6v. The quan- 
tity rjn represents the energy growth of a mode over one 
period - its inverse is equal to the quality factor (5„ com- 
monly associated with resonant electronic devices. The 
quantity /„ is the moment of inertia of the mode: 



|i5r„|^ dm 



(16) 



(e.g. Glasner & Buchler 1993). A positive work-integrand 
means that locally internal energy is converted into kinetic 
(pulsational) energy. 

In Fig. 12 we exhibit the work-integrands for the lowest 
14 modes as a function of z (The areas under the curves 
are the relative growth-rates). It is apparent in Fig. 12 that 
the egregious growth-rate of the strange mode is the result 
of two effects. First, the inner pulsation amplitude that 
we have shown to be characteristic of the strange mode 
causes both Ap and Ap and hence the integral (eq. 15) to 
be much smaller for the strange mode. Thus the surface 
region plays a relatively larger role in the growth-rate. 
Second, the smallness of the integral is compensated for by 
the considerable amplification that comes about from the 
small moment of inertia J„ in the denominator of Eq. 15. 
The resulting strange growth-rate is therefore largely the 
result of the contribution from the outer region. In this 
example, the outer region happens to be driving, this is 
the source of the positive growth rate of the strange mode. 

Since for the strange modes almost all of the driving 
comes from the HPIR the stability of these modes is there- 
fore more sensitive to the physical and numerical resolu- 
tion of the outer layers. However, as we have seen, their 
existence and nature are not sensitive at all. 



More light can be shed on the behavior of the growth- 
rates by examining a non-adiabatic extension of our toy 
model, which we consider next. 



6. Nonadiabatic Toy Potential Barrier Problem 

The nonadiabatic problem is of course much more com- 
plicated, but it is possible to extend the potential barrier 
problem of §2 to mimic nonadiabatic pulsations by gener- 
alizing Eq. 9: 



-y" + V k5{z - Zo)y + iUy 



Q'^y 



(17) 



where k is the (complex) wave-vector k = ^ (uj^ — iU\ 
and (i is the complex frequency. The potential \J mim- 
ics an overall dissipation. In the adiabatic case, V was a 
real constant of value Vo- Here we extend the potential 
to include an imaginary component V = Kj -|- iAVo- The 
choice of A14 represents a crude solution of the nonadia- 
batic pulsation equations: 



a;^(5r = Gi^r -h G2m 



(18a) 



iCoTbs = Ki^r Yi^Tis , (186) 

in which we assume that (G2(itI3 — Y^'2)~^Y^\)8r can be 
replaced with a complex term ;B(5r, proportional to br and 
thus to y. It has been recognized that the shape of the 
entropy fluctuations is independent of the mode (Pesnell 
& Buchler 1986) and is well represented by a sharp peak 
in the ionization region. As a result we superimpose the 
complex AVq onto the adiabatic delta function potential. 

The inner boundary condition is still y(0) = 0, but 
the outer one can now be more general (e.g. Morse & Fes- 
hbach, 1953, Eq. 11.1.21 ff), viz. 5p = Rii + Xu which 
translates here into y'{l) — Zy(l), where Z is a complex 
surface impedance. Z = corresponds to a free surface 
as in the adiabatic case (perfect reflection). ImZ is the 
inertia presented to the stellar surface by an overlying 
atmosphere, and ReZ corresponds to acoustic radiation 
losses. 

For highly nonadiabatic stars the adiabatic potential 
barrier that is contained in Gi becomes less important 
relative to the seond term, and it is the S function from 
the peaking of ds that dominates the problem. We note 
in passing that Zalewski (1992) decided to include an in- 
homogeneous term AS{z — Zq) instead of a potential, and 
took Z to be a real constant vq (whose value he did not 
specify in the paper) . His simple model does not work for 
Ccpheids although it matches some, but not all the prop- 
erties of the spectrum of the highly nonadiabatic PAGB 
models. 

With the appropriate matching conditions y{z^) = 
y{z:^) and y'{z+) = y'{z^) + kVy{zo), one obtains the 
eigenvalue equation for the complex wave- vector k 



cosfc —Zsmk 

= —Vo sin^i;o(cosfc(l — Zo) —Z sinfc(l — ^o)) 



(19) 
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Fig. 13. Period ratios (top) and relative growth-rates (bottom) 
as a function of outer well width 1-Zo for the lowest 9 modes 
of the complex barrier problem with V = 10 — 2i, U= 0.02. 



The solution is given by 

y{z) = a sin kz z < Zq 

= cos k{l — z) — Z smk{l — z) Zo < z 



(20) 



where the inner amplitude a{k) — y'{0)/ky{\) is now 
a{k) = (cos fc(l — Zo) — 2 sin k{l — Zo))/ sin kzo (21) 

Again, we examine a sequence of toy models, but this 
time for a relatively narrow range of barrier locations, Zo- 
We have chosen the parameters V = 10 — 2i and U = 0.02. 
We set = consistent with w 0, as in the hydrocode. 
The period ratios for this sequence are shown in Fig. 13 
(top) and it is easy to follow the strange property from 
overtone 9 at the bottom left corner to overtone 1 at the 
top right. In Fig. 13 (bottom) we show the corresponding 
growth-rates of the lowest modes, as obtained with Eq. 19. 

Considering the crudeness of the toy model's S func- 
tion potential, the properties of the growth-rates mimic 
those of the Cepheid sequence remarkably well. Here, as 
in Fig. 11, we see the increasingly rapid switchover of the 
unstable excursions and even, at bottom right, the weaker 
excursions of the second and third order strange modes. In 
this toy model, the effect of second and third order strange 



modes are also visible in the period ratios (Fig. 13 top) as 
it was in the Cepheid sequence (Fig. 11). 

The behavior of our eigenvalues can actually be un- 
derstood analytically via perturbation theory when the 
nonadiabaticity is small. The correction Aw^ (where a;^ = 
Lu'^ + Aw^ to first order) to the eigenvalue due to the com- 
plex potential U and the small imaginary contribution to 
V is given in terms of unperturbed (adiabatic) quantities 



iU + iAVoky^Zo)/! 

iU + iAVoa{kf sir? kzo/l 



(22) 



where I = J\y\^dz. The correction is seen to be at a 
maximum for the strange mode for which T is at a mini- 
mum. In contrast, it is at a minimum for the neighboring 
matched mode for which y{zo) '^0. Ths explains the rapid 
variation of the growth-rates in the vicinity of the strange 
modes Fig. 13 confirms both of these predictions. The po- 
tential U by construction provides an overall shift of Imw 
to negative values. 

7. The Cepheid Horn 

The reduction of the adiabatic pulsation problem to a 
Schrodinger like equation is very well known in the study 
of wind instruments (Morse & Ingard 1968, Benade 1977). 
But it is common to focus on pressure perturbations, 
though it is straightforward to show the equivalence of the 
following two descriptions of linear sound waves in a horn, 
the first for the velocity perturbation u (or cquivalently for 
5r), and the second, well known Bernoulli- Webster equa- 
tion, for the pressure 5p 



}_d_ ( P^d_)o ^ 
pdx\S dx^ ' 

S dx \p dx 



2 

OJ u . 



(23) 
(24) 



In the case of an acoustic horn it is not usually necessary 
to consider gradients in the sound speed Cg and density p, 
and it is the shape of the spatially varying cross-sectional 
area S{z) that determines the potential V{z). 

Eq. (23) can be transformed into a Schrodinger equa- 
tion of the form (Eq. 7) in the same way as Eq. (1), in 
which the potential is 



V^{z) = iS/p)-^{S/p)-^ 



(25) 



Mutatis mutandis, the potential for Eq. (12) is shown to 
be 



Vp{z) = iS/p)-^ 



(26) 



The potential Vp{z) is the so called iiorn function (Morse 
& Ingard 1968, Benade 1977). (We note in passing that a 
standard reference to the Webster-Bernoulli equation, viz. 
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Eq. 1 in Eisner (1967) incorrectly uses Eq. (24) instead of 
Eq. (23) as the equation for the displacement.) 

In the case of the horn it is the curvature asso- 
ciated with the flaring that can give rise to potential 
barriers. Reflection at these barriers (often phrased in 
terms of acoustic impedance mismatch) produces inter- 
nally trapped modes, which are the discrete tones of the 
horn or wind instrument in question. In contrast, loud- 
speaker horns are flared exponentially so as to avoid po- 
tential barriers. 

Can wc design a wind instrument with a flare function 
TZ{z) = \/ S{z) chosen so as to have the same acoustic 
spectrum as our variable star? In the Cepheid, as we have 
seen, the potential is produced by the variation of p, Cg, 
the spherical geometry, and by gravity. Can we work back- 
wards from the potential of Eq. (8) to flnd the shape 'R-{z) 
of horn that would reproduce a Cepheid's tones? 
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Fig. 14. Cross-section of the horn corresponding to Cepheid 
pulsations, shown over the outer potion of the pulsating enve- 
lope (The entire model envelope extends to « = 8); the kink at 
z = 0.75 corresponds to the potential barrier. 

A closer inspection of the stellar problem shows imme- 
diately a big difference with a horn, aside from the strong 
spatial gradients of density and sound speed: The poten- 
tial V in Eq. (8) contains a gravitational term that cannot 
be absorbed away into the term that has the form of a 
horn function. The presence of the gravitational potential 
term thus makes the inversion from the potential to the 
flare function impossible unless a Cowling approximation 
is made. 

With the Cowling approximation we can construct an 
approximate equivalent acoustic horn - with the horn part 
being of primary importance everywhere except in the in- 
nermost regions of the star. In Fig. 14 we show a slice 
through the horn that would produce the same horn func- 
tion as the Cepheid model (with Te//=5520 K). The modes 
produced by this horn of course have somewhat diflferent 
periods from the original Cepheid because of the Cowling 



approximation. It is indeed well known that this approxi- 
mation is not very accurate for low radial modes, and the 
reason has already been mentioned: The nodal structure 
is determined by the whole envelope, as the transforma- 
tion to z has shown particularly clearly. Nevertheless, it is 
amusing to see that in the region where the horn function 
dominates the potential, a constriction in the bell of the 
horn arises because of the HPIR. 

Because the strange modes are trapped surface modes 
they should feel very little the gravitational potential. 
They can therefore well described with the horn poten- 
tial only. Interestingly the radial displacement eigenvec- 
tors have an appreciable amplitude only the outer, open 
part of the flare. 



8. Nonlinear Pulsations 

The strange modes are seen to be unstable to the left of 
the blue edge of the fundamental mode and the first over- 
tone - in fact the strange mode can be the only unstable 
vibrational mode there (Fig. 11). It is of course interest- 
ing to investigate what the pulsations of this star look like 
and why no observational evidence has been reported for 
variability in this region of the HR diagram. 

We have seen in Fig. 8 that the region outside the 
HPIR dominates the radial eigenvectors and thus the pul- 
sation. A priori we have no idea of what the amplitude 
of pulsation might be. We have therefore computed the 
nonlinear (limit cycle) pulsations of several models from 
our sequence, that are located to the left of the Cepheid 
instability strip (and are therefore not Cepheids stricto 
sensu). 

For all these models the limit cycle pulsations have a 
very small amplitude, typically 0.1-1 km/s photospheric 
velocity and magnitude variations of milli-mags. If such 
'strange Cepheids' indeed exist they are very hard to de- 
tect. Aikawa & Sreenivasan (1995) have found similarly 
small amplitudes for yellow supergiant and intermediate 
mass stellar models pulsating in a strange mode. 

The effects of the strange mode might also be felt di- 
rectly in the regime of 'normal' Cepheids. Two simple sce- 
nari are possible here in theory: (a) the strange mode is 
linearly unstable together with the fundamental (or with 
the first overtone); if the nonlinear coupling coefficients 
happens to have suitable values, then beat pulsations can 
occur (The conditions that must be satisfied are described 
in great detail in Buchler & Kovacs 1986); (b) the strange 
mode is resonant through a higher order resonance, say 
4:1 or 5:1, cf. Fig. 11, top). If the strange mode is further 
linearly unstable or just marginally stable and the non- 
linear coupling coefficients are right, then periodic, fre- 
quency locked pulsations can occur (irregular amplitude 
modulations are also possible in principle, but less likely). 
Finally, the strange modes could have an indirect effect 
in that they could destabilize the fundamental (or first 
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overtone) limit cycles in a regime where the first overtone 
(or fundamental) limit cycles are also unstable. The result 
would then be a usual beat Cepheid, i.e. constant ampli- 
tude pulsations in both the fundamental and first overtone 
(Buchler 1997). 

For the strange modes to be able to have any of the 
described effects depends on the strength of the nonlinear 
coupling coefficients. Whether the small overlap between 
the normal and the strange modes allows a strong enough 
coupling is unknown at present. Hydrodynamic computa- 
tions are in progress to explore these phenomena. 



9. Conclusions 

We have shown that the problem of adiabatic radial pul- 
sations can be very instructively recast in terms of a 
Schrodinger equation - this without any approximations. 
The sharp variations in sound speed through the hydro- 
gen partial ionization region create a very high and very 
narrow potential barrier, separating the pulsating enve- 
lope into two regions. The appearance of strange modes is 
a consequence of the possibility of mode-trapping in the 
surface regions. 

While this Schrodinger equation formulation is stricty 
correct only in the adiabatic limit, the same phenomenon 
of trapping persists when nonadiabatic effects and acous- 
tic radiative losses are taken into account. 

It is clear from our study that the strange modes are 
thus of acoustic origin, rather than arising as a coalescence 
of two thermal or secular modes. 

The phenomenon of trapping has been illustrated with 
a simple analytical toy problem that exhibits all the im- 
portant features associated with the strange modes, such 
as avoided level crossings. This model can also be extended 
to explain the properties of the nonadiabatic spectrum. 
The resemblance with the behavior of the periods and 
growth-rates of a Cepheid model sequence is striking. 

There is of course no reason why strange modes should 
not occur in the nonradial spectrum of low I p modes, 
because these modes feel essentially the same potential as 
the radial modes (although a Schrodinger formulation is 
no longer possible without approximation because of the 
higher order of the nonradial equation) . The level crossings 
are therefore distinct from those that occur when p and g 
modes interact (Unno et al. 1989). 

Finally, we have examined the nonlinear pulsations as- 
sociated with the strange modes and have concluded that 
they have very small amplitudes, especially in luminosity. 
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